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Abstract 

We analyze the transverse structure of the Gouy phase shift in light beams carrying orbital 
angular momentum and show that the Gouy radius re characterizing the transverse structure 



grows as y^2p + |^| + 1 with the nodal number p and photon angular momentum number i. The 
Gouy radius is shown to be closely related to the root-mean-square radius of the beam, and the 
divergence of the radius away from the focal plane is determined. Finally, we analyze the rotation 
of the Poynting vector in the context of the Gouy radius. 
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I. INTRODUCTION 



Light beams carrying orbital angular momentum (OAM) jl-|3|, which is distinct from spin 
(or polarization) angular momentum, have become important subjects for theoretical and 
experimental investigation over the past two decades, and have been studied in the context 
of applications as diverse as angular momentum transfer to trapped particles , quantum 
communication ^6 8 ] , and sensing Q. The OAM is directly related to the azimuthal phase 
of the mode, independent of the axis about which the phase angle is defined 

0L 

Photons 

exhibiting orbital angular momentum are also referred to as twisted photons 



collimated coherent beam of such photons is called an optical vortex 



n 



ii 



ll2j. A 



or a helical laser beam 



13| , a twisted beam [14( , 
15| . The interaction of beams carrying OAM with matter has become 



a subject of considerable interest 



ifHiai 



In the paraxial approximation, cylindrically symmetric light beams carrying OAM may 
be described by Laguerre-Gaussian (LG) modes, which are a natural generalization of the 
pure Gaussian beam without OAM. Laser beams with a pure Gaussian transverse profile 
are characterized by a Gouy phase of the form 

exp[-z0 o ,o] = exp[-itan _1 (2;/^)] , 

where z is along the direction of propagation, zr is the Rayleigh range (to be defined below), 
and the focal point or beam waist is taken to be at z = 0. It is evident that as the beam 
propagates from — oo to +oo it undergoes a phase shift of 7r relative to a plane wave [20 ]. 
A modified Gouy phase appears for an LG beam with winding number £ and radial node 
number p > 0, 



exp[— i9 P) g\ = exp 



-i(2p+ \l\ + l)tan- 1 (j-j 



Thus the corresponding cumulative Gouy phase change of an LG beam is ir(2p + \£\ + 1). 

Historically the Gouy effect has been viewed primarily - if not exclusively - as a longi- 
tudinal phenomenon, with the beam accumulating a phase change of 7r, or more gene rally 



, , , m ,,,,,,,, .fly 
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22]. 



However, there is also a transverse Gouy effect that has received scant attention. 

We will show that the Gouy phase shift in an LG beam has a transverse radial dependence. 
The Gouy effect is maximal along the axis of propagation, conventionally taken to be the 
z axis, and decreases radially outwards. At a specific value of the radial coordinate r, the 
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Gouy phase change vanishes, and it becomes negative for higher values of r. We shall call 
this critical value of r the Gouy radius and label it as ra- 
in what follows we shall establish that: (1) In the focal plane of a pure Gaussian beam, 
the Gouy radius ra equals the beam waist w(0). (2) In the focal plane of an LG beam, 
the Gouy radius r G becomes w(0)^2p + \£\ + 1, and is proportional to the rms radius of 
the beam, which also grows as ^2p + \£\ + 1. (3) The Gouy radius tq increases with z, 
becoming divergent at the Rayleigh range Zr, even while the width w(z) and the rms radius 
r rms (z) remain finite. (4) The Gouy radius serves as a natural physical boundary for the 
photon beam, delimiting the "allowed" zone from a "forbidden" zone. (5) The Gouy radius 
is helpful in understanding the rotation of the Poynting vector in a twisted photon beam. 



II. LAGUERRE GAUSSIAN BEAMS 



A beam of photons propagating in the z direction prepared with linear momentum fc , 
ar momentum £, and radial node number p may be represented by the vector 



orbital angu 



potential [23 1 



A = xw(r, z,0) < 



ik()Z 



(1) 



LG Pi i(r ,z,' 



2p\ 



where x is a unit vector in the x direction and 

u(r,z,(j)) 

x exp[i£(f)] exp 
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25j: 



Tt(p + \£\)\ w(z 
ik r 2 z 



L2(^ + 4)j 



cxp 



V2r 
w(z) 

i(2p+\i 



exp 



w 2 (z) 



L Wl 2r2 

p yw 2 {z) 



1) tan -1 — 



(2) 



where Lp are the generalized Laguerre polynomials, and the Rayleigh range zr will be 
defined below. For a pure Gaussian beam (p — I — 0), the radial dependence of the 
amplitude is given simply by 



u(r, z, <p) ~ exp 



w 2 (z) 



(3) 



and the beam width or radius of such a beam is conventionally defined as the radial distance 
at which the vector potential falls to 1/e of its value on the central axis r = 0. Thus by 
definition the width or radius of a Gaussian beam is w(z). 

A general LG beam, however, has an additional radial structure, including p radial nodes. 
Yet the parameter w(z) that appears in Eq. (T5]) is still conventionally referred to as the width 




G Negative Gouy effect 
Zero Gouy effect 
Positive Gouy effect 

p Maximum Gouy effect 

Positive Gouy effect 
Zero Gouy effect 
H Negative Gouy effect 



OA = OB = Gouy Radius 



FIG. 1: Section of a Gaussian Beam. 
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281 ] . But this is misleading for nonzero i or p, as 



or radius of the LG beam 
we shall demonstrate below. 

At the beam waist of a Gaussian beam, where the radius is w(0), the cross sectional area 
is A = 7tw 2 (0). A is also the area of the aperture from which the Gaussian beam emerges. 
A quantity of importance for laser beams is the Rayleigh range Zr = A/X, where A is the 
wavelength of the beam (Ref. 29J, p. 714). So an important relationship exists between the 
beam waist radius w(0), the Rayleigh range zr and the wavelength A of a Gaussian beam 



29): 



w 2 (0) 



Xz R 



Ti 



(4) 



Equation H] may be used as an alternative definition of the Rayleigh range for a Gaussian 
beam. The Rayleigh range sets the scale on which the Gaussian beam spreads 29], 



w(z) = w(0) 



z 2 + zl 

4 



1/2 



(5) 



III. TRANSVERSE GOUY EFFECT FOR GAUSSIAN AND LG BEAMS 

Because of the divergence of the Gaussian beam as it propagates away from the beam 
waist at z = 0, the wavefront is convex rather than planar, as illustrated in Fig. [TJ So the 
longitudinal wave vector component k z is smallest along the z axis (where the Gouy effect 
is strongest and the phase advances most slowly) and increases with the distance from this 
axis. Of course, for a perfectly collimated beam with zr = oo, the wave front would be 
planar, and k z would be everywhere a constant equal to ko- But for a Gaussian beam k z is 
not constant across the transverse profile. One might expect that ko would equal the value 
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of k z at the center of the beam, on the z axis. But this is not the case. Because of the Gouy 
effect, the value of k z on the z axis is necessarily smaller than the wave number ko which 
would correspond to a plane wave front. 

In what follows we shall examine the Gaussian beam at its waist, at z — 0, and show 
that it is precisely on the circumference of the beam waist, at r = w(0), that k z = ko, and 
that k z increases monotonically with r. 

The normalized field pattern for a Gaussian beam is given by setting p = i = in the 
general expression of Eq. fl2]) for LG beams: 



u(r,z,<f>) 



2 exp[ik z — ii^(z)] 



exp 



+ 



ik n r 2 z 



w 2 (z) 2(z 2 + z 2 R )_ 



(6) 



7T w(z) 

where ip{z) = tas\~ x {z/ zp). Now k z can be evaluated as a function of the radial distance r at 
the beam waist z = 0, by obtaining the derivative —idu/dz. To first order in z, ip(z) — > z/zr, 
and from Eq. (jSJ) w(z) — > w(0), so 



kJr 



z=0 



It is evident that k z = ko when 



k + 



2z r z r X 



1 k r 2 



2z\ 



(7) 



ko vr 

where A = 2-7r//co is the photon wavelength. Clearly, the Gouy radius at z = is rc(0) = 
w(0). So we have established the important result that for a Gaussian wave at its waist or 
focus, the Gouy radius is identical with the beam waist radius w(0). 
Similarly for the general LG beam (Eq. (J2J), we have 

kor 2 -2(2p+\£\ + l)z R 



k z (r) z=0 = k + 
This yields the Gouy radius 

r£(0) = (2p+ |*|+1) 



24 



(8) 



Xzr 



71 



{2 P + \e\ + i)w 2 {o) 



(9) 



IV. RADIAL INTENSITY PROFILE OF TWISTED BEAMS 



The amplitude of the LG mode is given by 

\u(r,z,<j)) \ = 



2p\ 



tt(p + \e\y.w 2 (z 



. x W* e -*/2 L \t\( x ) 



where x = 2 f^ > . This function oscillates in a range x_ < x < x+, where 

x± = 2p + |£| + 1 ± v/(2p + \l\ + l) 2 - Pi - l) 2 (10) 

and decays exponentially outside this range [yj]. We then have 

C°° m 2 (z) 
(r 2 )= / 27ir i dr\u(r,z, ( p)\ 2 = ^-t(2 P +\£\ + l), (11) 
Jo * 

and 



r rms (0)=^ 2p + jf l + 1 . (12) 

Thus for a pure Gaussian wave, with p = i = 0, r rms = w(z)/y/2. Notably, we find that at 
the beam waist z — 0, the rms beam radius is l/v2 times the Gouy radius computed in 
Eq. (jHJ), for all values of p and |£|. 

Clearly, r rms is not the same as the radial distance of maximum intensity r max which has 



the value (for p = 0) [31 ] 



r m axW=^)\/y- (13) 

The radius of maximum intensity becomes zero for a pure Gaussian wave (p = £ = 0), but 
the expected or rms value of r does not vanish for any realistic distribution. For large values 
of £ at fixed p, r max approaches r rms . 

V. BEAM WIDTH AWAY FROM THE WAIST 

We shall now tackle the more general problem of finding the beam width of a twisted 
photon beam for arbitrary values of z, not just at the beam waist or focus. A note of 
caution is in order, because Eq. (J2J) is a solution to the Maxwell equations only in the 

beam the equation is physically meaningful 
only for values of z sufficiently small in relation to the Rayleigh range. We will write Eq. (j2J) 
in compact form as 

u = F(r, z) exp iG(r, z, 0) , (14) 
where F and G are real. So using Eq. ([1]) we may write 

i du 

k z = k — (15) 

u oz 
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As before, we will obtain the wave vector or momentum component k z , this time as a function 
of both r and z: 

.du .IdF dG 
- l ^z = - l F^ U + ^ U - (16) 

The photon momentum is the sum of an imaginary and a real term. The former corresponds 



to an amplitude change, and the latter to a phase change [2l|]. We are interested in the phase 
change, so 



7„ , dG fc r 2 (4-^)-2^(4 + ^)(2p+|£| + l) 
kz ~ ko + -dz~- ko + 2(4 + z^2 • ( 17 ) 



Using the definition of the Gouy radius tg{z) as the value of r for which this change in phase 
becomes zero, we have 

rUz) = ^(2p+\e\ + 1)^5 = w\0)(2p+ \e\ + 1)^5- (18) 

71 \ Z R z ) \ Z R z ) 

This expression reduces to Eq. (E]) as z — > 0. The Gouy radius also diverges as z — > Zr. But 
as stated earlier, this latter case is not physical, as the paraxial approximation of Eq. (J2]) 
breaks down for large values of z. For \z\ Czjjwe may simplify the result as 

rU z )=w 2 (0)(2p+\i\ + + . (19) 

VI. TOTAL LINEAR MOMENTUM OF A LAGUERRE GAUSSIAN BEAM 

The linear momentum of a Gaussian wave photon differs from the corresponding momen- 
t U ,„ of a p.a„e wave photon because of the Gouy phase shift fl, as is evideu. tau Bq. ® . 
We will now calculate the linear momentum of a Laguerre Gaussian photon. 

We first find the average momentum density integrated over the beam waist, i.e., over 
the z = focal plane. From Eq. (jSJ) we obtain 



k , 2v 2p+ \£\ + 1 

; z/2=o ''0 • 

and using Eq. (II ID . we have 



(k z } z=0 = k + -\(r 2 } - -^-^ , (20) 

IZ R Zr 



{k ^ s=kls .^±mi. (21) 

For a pure Gaussian beam, this reduces to (k z ) z=0 = k — 1/2zr. For a plane wave, for which 
Zr = 00, the photon momentum reduces to k . 



To obtain the average momentum density at arbitrary z, we use Eq. (fITj) and average 
over a plane normal to the z axis (the propagation direction): 

{K) = fa + _ ^ + 14 + = fa - (^±1) . (22) 

2(z^ + z 2 ) 2 z 2 + z|, 2zh 

We see that the result is ^-independent, expressing the conservation of photon momentum 
along the direction of propagation. 

VII. PHYSICAL SIGNIFICANCE OF THE GOUY RADIUS 

The (longitudinal) Gouy effect has been exp lained as a result of the transverse confine- 



ment or restriction of the photon beam 
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30|. This restriction leads to the generation of 



transverse photon momentum or wave vector components, so that 30] 



{k.) = ko- { ®T { ® , (23) 
where it is to be noted that (k z ) 2 7^ (kf). Comparing with Eq. ( 12"2"]1 we obtain 

<*» + tfl = 2p { V { o\r 1 • (24) 

where the average is performed over a plane of constant z, or, using the cylindrical symmetry 
of the beam, 

{k x )-{k y }- 2(w(Q))2 • (25) 

These equations express the position-momentum uncertainty in the transverse direction. 

Now if we consider averaging over only in Eq. ( |23l) . at fixed z and r, we see that 
within the Gouy radius the transverse wave numbers k x and k y are real, but beyond this 
radius k x and k y become imaginary. Imaginary wave numbers imply amplitude damping. 
The situation is analogous to a particle confined by a finite potential wall. There is a finite 
probability of penetration within the potential wall, and the probability drops off rapidly 
within this "forbidden" zone. Thus the Gouy radius provides a more satisfactory "physical" 
boundary for the photon beam than the function w(z) or the related radius of maximum 
intensity or the rms radius. 

We shall now offer a correction to an important result that was derived for the rotation of 



the Poynting Vector in the course of the propagation of a Laguerre Gaussian beam 31] . The 



authors have obtained an expression for the winding of the Poynting Vector around the axis 
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of propagation of the photon beam. Taking the Poynting Vector to be concentrated at the 
radius of maximum field they showed that this Poynting Vector undergoes a rotation of 71/ 2 
when propagating from the beam waist to the far field, and that the rotation is independent 
of the value of £. But a more accurate calculation shows that the Poynting Vector does 
indeed circulate as a corkscrew and that the rotation is dependent on the value of £. 

The Poynting vector S of an LG beam is obtained by taking the gradient of u in Eq. fl2]): 



S 



c r, z) 



zr 



z 2 + z 2 



:f + 



1 



T {Zr — z ) 

2{z 2 R + z 2 ) 2 



(2p+\£\ + l)z R 
(4 + * 2 )ko 



(26) 



where H is a function of r and z. (The expression given in Ref. 3l| is missing a square and 
a minus sign.) In Ref. 3JJ, the second and third terms inside the parentheses are assumed 
to be small and neglected. But for r of order of the beam width, r 2 ~ z R /k Q (Eq. (j3J)), these 
terms contribute to the Poynting vector magnitude \S\ at the same order as the term. 



to 



What is clear is that dropping these two terms is equivalent to setting their sum equa 
zero, which is valid but only at the Gouy radius ra- So in the expression given in Ref. [31 
for the rate of azimuthal rotation, 



dz 



k r 2 



(27) 

t is therefore problematic to 



r is not an arbitrary radius but precisely the Gouy radius Tq. 
assign to r the value of the radius of maximum field intensity [3jJ or any radius other than 
the Gouy radius. 

Substituting for tq from Eq. ( Fl8|) . we obtain 

1*1 (4 - z 2 ) 



dip 
dz 



2z R (2p+\£\ + l) (z 2 R + z 2 



(28) 



and upon integrating, we have 



(2p+\£\ + l) 



tan 1 ( — 



z 



2z R 



(29) 



Thus far, we have not made any approximations. Eq. ( 129]) follows exactly from Eq. (J2J). But 
Eq. 02]) itself is an exact solution to Maxwell's equations only in the paraxial approximation, 
z z R . Eq. (128]) shows that at small distances from the focus, z z R , the rotation rate is 



dy _ \e\ 

dz 2z R {2p + \£\ + 1) 



(30) 



to be compared with the rate dip/dz = l/z R obtained in Ref. [311 ]. 
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VIII. CONCLUSIONS 



We have analyzed the transverse structure of the Gouy effect in a Laguerre-Gaussian 
twisted photon beam, with an emphasis on the Gouy radius r G , defined as the critical 
radius where the Gouy phase shift vanishes. We have shown that in the focal plane of the 
LG beam, the Gouy radius is proportional to the rms radius of the beam, which in turn 
scales as ^/2p + \£\ + 1 with the angular momentum £ and nodal quantum number p of the 
beam. However, away from the focal plane, vq grows faster than the rms radius, and diverges 
at the Rayleigh range zr, even while the beam width and the rms radius remain finite. The 
Gouy radius is shown to delineate a boundary between the "allowed" and "forbidden" zones 
of the photon beam, and is conducive to understanding the rotation of the Poynting vector 
in a twisted beam. 
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